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ABSTRACT 

Aims. We study two dimensional Fabry-Perot interferometric observations of the nearby face-on late-type spiral galaxy, NGC 628, in 
order to analyse the ionized gas component of the interstellar medium. Covering the galaxy out to a radius larger than 12 kpc, and 
with a spatial sampling of r.'6, we aim to investigate the large-scale dynamics as well as feedback from individual Hii regions into 
their surrounding medium. 

Methods. The observed Ha emission distribution and kinematics are compared with auxiliary data from molecular and atomic gas 
observations, which display many similarities. We decompose the observed line-of-sight velocities into rotational and higher-order 
harmonic components, and study the role of gravitational perturbations along with that of external triggers which can disturb the 
kinematics and morphology of NGC 628. We calculate radial profiles of the emission-line velocity dispersion which we use to study 
the role of feedback from individual Hii regions. 

Results. We verify the presence of an inner rapidly rotating disc-like component in NGC 628, which we interpret as caused by slow 
secular evolution of the large-scale spiral arms and oval structure. In combination with auxiliary data, we find indication for that gas 
is falling in from the outer parts towards the central regions, where a nuclear ring has formed at the location of the inner Lindblad 
resonance radius of an an in = 2 perturbation. Complementary continuum subtracted narrow band images in Ha have been used to 
identify 376 Hii regions with calibrated luminosities*. The mean velocity dispersion for the ionized gas (even when excluding pixels 
belonging to H ii regions) is almost constant out to 12 kpc, although it varies from 14 to 20 km s"', with a steady decline in the outer 
parts. 

Conclusions. We have found kinematic signatures of radial motions caused by an m = 2 perturbation. Such a perturbation may well 
be responsible for the inflow of material forming the nuclear ring and the inner rapidly rotating disc-like structure. The latter, in turn, 
could help build a pseudo-bulge in NGC 628. The current paper demonstrates a number of tools that we have developed for building 
a solid frame work for studying the evolution of structure in spiral galaxies using two dimensional kinematic observations. 

Key words. Galaxies: spiral ~ Galaxies: structure - Galaxies: evolution ~ Galaxies: kinematics and dynamics - Galaxies: ISM - 
Galaxies: individual; NGC 628 



1. Introduction 

Constraining the three dimensional shape of galaxies is one 
of the most important elements in validating galaxy formation 
and evolution theories. While images provide projected mor- 
phological information, complementary kinematical data is cru- 
cial to constrain the shapes and the physical processes that 
govern them. The vertical components of velocity and veloc- 
ity dispersion (cf) are perhaps the most important parameters 
required to determine t he flattening of the various disc lay- 
ers (e.g.. lBottemall2062h. the formation and evolution of bars, 
bulges, and discs (e.g.. lCombeslll99ltlvan deii Boschll 19981). the 
intrinsic shape of the dark matter component ( Zasov et al.lll99ll: 
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* The complete version of the Hii region catalogue. Table 1 in the 
online section of this manuscript, is only available in electronic form at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-l-A/| 



iBecquaert & Conibe s|ll997| ), and the nature and extent of disc- 
halo interactions ^Norman & Ikeuchiiri 988). 

If the gaseous layers in the disc are isothermal, both atomic 
and molecular gas should have similar heights. Accordingly, the 
atomic and molecular components can be considered as two dif- 
ferent phases of the gas coexisting as the same dynamical com- 
ponent, and the vertical oscillations will be of the same ampli- 
tude. It is possible that the density of clouds and/or the radiation 
field at high altitude is not enough to excite the CO molecule, 
which means that the limit for observing CO will not be coin- 
ciding with the limit for neutral hydrogen itself. There is a sharp 
boundary where the apparent molecular CO fraction falls to zero, 
while we expect a smoother profile for a unique dynamical gas 
component. However, due to a lack of kinematic observations of 
H 11 regions in galaxies, previous studies of the molecular com- 
ponent have still not been able to verify whether we can observe 
CO formation and destruction, and therefore put constrains on 
the velocity and velocity dispersion of the different gaseous com- 
ponents. 
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Alternatively, the gas is not isothermal and several layers of 
gas with different compositions or different thicknesses coexist. 
The mixing time scales, through atomic or molecular collisions, 
of two different layers obey the second law of thermodynam- 
ics. For a diffuse and homogeneous gas the mixing time scale 
is of the order of 10"* years. In more realistic cas es, the colh - 
sional time scale is of the order of 10 years (e.g., Cowid ll980l) . 
This value is comparable to the dynamical time, but is much 
larger than the chemical mixing time scale due to the phase 
change of the gas, estim ated to be of the order of 10^ years 
jLanger & Graedelll989l) . Implicitly, the basic constituents of 
the interstellar medium, i.e., molecular or atomic clouds, cannot 
be considered as a fluid in equilibrium. If collisions redistribute 
the kinetic energy, the intrinsic velocity dispersion would be in- 
versely proportional to the cloud mass, whereas in the case of su- 
pernova feedback, the resulting gas velocity dis persion is in th e 
range of 4-7 km s"' throughout the stellar disc (lBottemall2003l) . 
T hus, the gas can be stirred up to higher dispersions as was noted 
bv lGerritsenI ( Il997h . Consequently, in the outer parts of the disc, 
where star formation is less, one might expect a decrease of ve- 
locity dispersion (Bland-Hawthorn e t al.ll 19971) . 

Repeatedly, this picture has been contradicted by observa- 
tions (e.g., Jimenez- Vicente et al. 1999). At large galactocen- 
tric radii, where streaming motions should be minimal, the cr 
for giant molecular clouds remains almost constant over sev- 
eral orders of magnitude in brightness. Thus, no equipartition 
of energy needs to take place and other mechanisms respon- 
sible for the heating are required. Relativel y small clouds can 
be he ated by star formation and supernovae jChieze & Lazarefl 
Il981h . whereas t he largest clouds coul d be heated by gravita- 
tional scattering (i Jog & Ostrikerl [r988l) . But collective effects 
such as gravitational instabilities forming structures such as spi- 
ral arms, have not yet been taken into account. 

In this context, it is crucial to compare the processes that 
govern the ionized gas with those that control the molecular and 
atomic gaseous components. Although the neutral and molecu- 
lar components are commonly used to derive large-scale kine- 
matic information, higher resolution ionized gas observations 
are necessary to better study the small-scale velocity structure 
due to recent star formation and corresponding feedback pro- 
cesses. Face-on galaxies are fav ourable targets to derive the lo- 
cal t urbulent velocity fields (e.g.. lShostak & van der Kruitll984t 
ICom bes & Becquaert 1997), since the line width can be at- 
tributed almost entirely to the velocity dispersion in the vertical 
direction. Previous studies have shown a remarkably constant 
(o-^ 7kms"') gaseous cr as a function of galactocentric radius. 
This implies that the vertical streaming motions at the arm cross- 
ing are not dominant kinematic features. In inclined systems, due 
to the augmented effect of non-circular mo tions in the plane of 
the disc, these studies are complicated (e.g.. lVogel et al.lll993[) . 

In this paper, we present two-dimensional kinematic maps 
of the Ha-emitting gas in the nearly face-on galaxy, NGC 628. 
We apply a detailed kinematic analysis to explore the role of 
large-scale versus local ionized gas distribution and dynamics. In 
Section|2] we present a brief review of relevant previous studies 
of this galaxy, and in Section [3] we describe our observations. 
In Section ID we present our kinematic analysis tools as well as 
the results, and finally in Section|5]we discuss our results in the 
context of secular evolution of discs in late-type spiral galaxies. 
Finally, in Section|6l we conclude the implications of our results. 



2. NGC 628 



NGC 628 (or M74) is an extensively studied isolated grand- 
design Sc spiral galaxy at a distance of 9.7 Mpc (implying 
1" ^ 47 pc), and with a B-band magnitude of Mb=-19.98 
(lTullvl[T988) . The inner disc of NGC 628 is relatively unper- 
turbed, with an inclination of 6.5° and a position angle (PA) 
of 25 ° (e.g., iKamphuis & Briggsl 119921 : ICombes & Becquaen 
Il997h . This galaxy has a Holmberg radius of 6' with a rela- 
tively regular neutral gas velocity field inside thi s radius (e.g., 
IWevers et al.|[l986.). From ultraviolet observations. lCornett et al.l 



19941) have not been able to confirm the presence of a bulge, but 



found instead that the nucleus of NG C 628 has the appearance 
of disc material. Natali et all d 19921) had earlier suggested that 
NGC 628 could be seen as an inner and outer disc, characterised 
by different stellar populations, with a transition radius located 
at about 1.5 Holmberg radii. Although the exact location of this 
transition can be questioned. ICornett et al.l (1 1994 1) confirmed this 
picture, finding that the star formation history of NGC 628 varies 
with galac tocentric distance. Near infrared images s how an oval 
distor tion dJames & Seigaii [19991) and a weak bar iLaine et al] 
l2002h . of radii 2 kiloparsec (kpc) and 100 pc, respectively. 

Many previous studies have focus ed on the neutral hydro- 
gen kinematics and distribution (e.g., IShostak & van der KruitI 
1984) . the nature, distribution, and composition of Hn re- 
gions (e.g.jBellev & Rovlll992h. and t he molecular gas compo- 
nent (e.g.. Combe s & Becquaertlll997' ). Neutral hydrogen stud- 
ies have revealed the presence of an elongated ring-like struc- 
ture at around 2 Holmberg radii (12') distance from the nu- 
cleus. This ring lies in a plane with ^ 15° in clination with re- 
spect to the plane of the bright inner disc (e.g.. lWakker & Adleil 
1995'). A strong kinematic PA twist of the H i velocity field in- 
dicates that the outer part of the disc, where the ring is ob- 
served, is warped. The apparent isolation of NGC 628 makes 
it difficult to pin down the origin of the warp. The asymmet- 
ric gas distribution in the outskirts of the disc, however, could 
indicate some recent disturbance, as such asymmetric distribu- 
tions cannot survive for a Hubble time. Alternatively, the detec- 
tion of two large high-velocity clouds, symmetrically located on 
opp osite sides of the disc, co uld explain the warped disc (e.g., 
Wa kker & van Woerden|[T99TI) . These cl ouds could be accreting 
onto t he outer parts, causing the warp (ILopez-Corredoira et al.l 
120021: iBeckman et al.ll2003.) . 

This galaxy was a lso observed with the two-dimensional 
spectrograph SAURON ( Bacon et aT]|2001h . Although the field of 
view of SAURON (33" x 41") is much smaller than those used in 
all the other studies mentioned above, this instrument delivers 
stellar as well as ionized gas kinematic information, simultane- 
ously. The stellar cr decreases from ^ 80 km s"' in the outer parts 
to ^ 30 km s"' in the central zones, indicating a <20" radi us dy- 
namically cold central region, maybe an inner disc (Gan da et alj 
2006 ). Outsi de this region, the SAURON stella r cr is similar to 
that found by Ivan der Kruit & Freemaiil ( |1984|) . The latter study 
shows values of 60 + 20 km s"' at about one luminosity scale 
length, and its evolution with radius is compatible with an ex- 
ponential decrease, with a radial scale length twice that of the 
density distribution. This is ex pected if the stel lar disc has a con- 
stant scale-height with radius (lBottemalll993l) . The ionized gas 
kinematics follows that of the stars. The H/? distribution is more 
extended than the [O iii] distribution, both suggesting a ring-like 
structure of about 15 -20" radius (Ganda e t_ani2006() . This sup- 
ports the results from lWakker & Adler, (il995l) . 
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3. Observations 

NGC 628 is a member of the large sample of galax ies that were 
presented in the Fabry-Perot kinematical survey of Daig le et alj 
(2006a). The galaxy was observed with the two dimensional 
spectrometer FaNTomM on November 1 8th 2003, on the 1 .6 meter 
telescope at Observatoire du mont Megantic in Quebec, Canada. 
FaNTomM includes a narrow-band (10-20 A) interference filter, 
a Fabry -Perot interferometer and a photon counting imaging 
device ( Hernan dez et al.ll2003h. The effective focal length over 
telescope diameter ratio calculated from the effective pixel size 
on the sky is 2.00. The pixel size is 1'.'6, with an effective field 
of view of about 13.7' x 13.7', thus the data cube delivered by 
FaNTomM has the spatial dimensions 512x512 pixels. 

Our observations of NGC 628 were carried out with a filter 
centred at Ac = 6598 A, with maximum transmission at Ac and at 
293 K of 73%, and with a Full Width Half Maximum (FWHM) 
of 18.2 A. We scanned 64 channels of 0.11 A each. This corre- 
sponds to a velocity step of 5.2 km s"' per channel. Each chan- 
nel exposure was for the duration of 2.33 minutes, adding to a 
total of 149 minutes effective exposure time. The Fabry-Perot 
interference order was tuned to p - 899 at the Ha rest frame 
wavelength, implying a free spectral range of333.5kms '. The 
Fabry-Perot finesse across the field of view was 23.6, yielding 
an effective spectral resolution of 21216. 

3.1. Basic Data Reduction 

Although the data for this ga laxy have been reduced and pre- 
sented bv Da igle et al.l (12006 a^). we re-reduced the data in a man- 
ner customised for the analysis presented in this paper. Stacking 
individual exposures of each channel results in three dimen- 
sional interferograms. In order to transform these raw interfer- 
ograms into wavelength-calibrated data cubes, we created two- 
dimensional phase maps using observations of the narrow Neon 
emission line at 6598.95 A. This line has been chosen as it 
falls close to the red-shifted Ha line, limiting the phase shift in 
the dielectric layer of the Fabry-Perot. We then applied a slight 
spectral smoothing to the galaxy spectra as well as the calibra- 
tion spectra by convolving with a Gaussian with a FWHM of 
10 km s"', followed by sky subtraction by assuming a uniform 
OH emission profile throughout the field. Thus, we extracted ap- 
proximately 20000 spectra, from regions with no Ha emission, 
evenly spread over the field, averaged all these spectra, and sub- 
tracted one representation of the OH emission from the entire 
data cube. The spectral smoothing is applied in order to improve 
the derivation of the FWHM, given that, in the low signal re- 
gions, non-smoothed data deliver overestimated widths. These 
part s of the reduction w ere carried out using the ADHOCw pack- 
age dAmram et al.lll998l) . 

An improved data reductio n and p reparation package has 
been presented by iDaigle et al.l ( l2006bl) . We tested the robust- 
ness of the parameters we have used by cross-checking our 
results with corresponding files obtained by the latter reduc- 
tion package, findin g no significant differences. In particular, 
iDaigle et al.l (l2006b l) introduce a two-dimensional adaptive bin- 
ning method. Although this results in increased S /N in the outer 
parts, we decided not to use the binning since it reduces the spa- 
tial resolution in too many regions throughout the observed field. 

3.2. Deriving the Ha Distribution and Kinematics 

We derived the kinematics by fitting one Gaussian profile to each 
individual emission line. The fits yield the continuum level, the 



monochromatic Ha intensity, the Une-of-sight velocity, and the 
FWHM for each pixel independently. We estimated the errors 
for the derived kinematics i n a rigorous way by using the boot- 
strap Monte-Carlo method (Press et aljfl992l) . Accordingly, we 
fitted the kinematics using Gaussian randomised emission line 
representations. The errors could be derived in a robust way by 
repeating this procedure 200 times. In regions with high S IN, we 
can derive the velocities to within 4 km s ' and the cr to within 
5 km s"'. These values increase to 7 km s"' and 8 km s"', re- 
spectively, when averaging over all pixels within our field. 

We applied a spatial smoothing by convolving each plane 
of the data cube with a Gaussian to reach a final effective spa- 
tial resolution of 3". Testing various final resolution values, we 
found that this values was a good compromise between spatial 
resolution and coverage, given the intrinsic resolution on the sky 
of the original data cube. We then calculated the instrumental 
dispersion using the weighted mean of 10000 calibration spec- 
tra, and found cri„str=(FWHM/2.35)=7.37 ± 0.14 km s"', which 
we subtract quadratically from the derived emission line cr in or- 
der to correct for the instrumental broadening. The low error on 
the derived instrumental dispersion demonstrates that this quan- 
tity stays virtually constant across the observed field. To remove 
any inadequate fits to the spectra, we then employed a number of 
selection criteria. We removed all pixels for which the amplitude 
over noise (mean standard deviation of the emission-free region 
of the spectra) was less than five, or the continuum level was 
less than zero, or the intrinsic velocity dispersion cr < 0. The 
final result from this cleaning procedure leaves us with 21952 
independent kinematic measurements within our field of view 
(see Fig.[T]). 

3.3. Extracting tlie Parameters of IndividuaiHii Regions 

Ha imaging of NGC 628 was obtained with the 1 m Jacobus 
Kapteyn Telescope on La Palma. The observations were taken on 
July 6 2003 with photometric sky and \'!2-V!3 seeing. A SITel 
2048 X 2048 CCD detector was used providing a projected pixel 
scale of O'.'33/pix and a field of view of ^11'. 

The galaxy was observed through a narrow band (63 A width) 
interference Ha filter for 2600 s (2 exposures of 900 s plus an 
exposure of 800 s) and through a Harris /^-filter for the contin- 
uum subtraction for 2x500 s. The width of the Ha filter allows 
for partial transmission of the nitrogen [Nii]/16548,6583 emis- 
sion lines. This contribution has not been subtracted from the 
measured fluxes, and is unlikely to amount to more than 20% of 
the total in any specific case. 

The data reduction was carried out using standard IRApQ 
tasks. They include over-scan subtraction and bias and flat field 
corrections. Afterwards the sky emission was removed from the 
individual images by subtracting a constant value obtained from 
measurements on areas of the image free of galaxy emission. 
The individual images for the Ha and /?-band filters were then 
aligned using the positions of non saturated field stars in both im- 
ages and combined separately. Then, the on and off -band images 
of the galaxy were astrometrically calibrated using the USN02 
catalogue coordinates for the foreground stars of the galaxy im- 
ages and the IRAF tasks ccmap and msctpeak. The accuracy of 
the astrometric calibration was ^0.34". 

Finally the continuum image was subtracted from the on- 
line (Ha-i-continuum) image. To estimate the scaling factor for 
the continuum subtraction we measured the flux of ^15 non- 



' IRAF is distributed by tlie NOAO, wfiicli is operated by AURA un- 
der cooperative agreement with tlie National Science Foundation. 
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Fig. 1. The distribution and kinematics of the Ha emission derived from our Gaussian fitting routine as labelled above the individual 
panels. The north-east direction is indicated at the bottom-right corner of the continuum panel. Also illustrated are the results from 
our various models for the observed Ha velocity field. All maps have been centred around the kinematic centre (0,0) derived in an 
iterative way described in Section |4!2] The exponential disc model with corresponding residual field (data - model) is presented in 
the middle row. The non-parametric purely rotating disc model, the harmonic reconstruction, and the residual (data - harm, rec.) are 
in the bottom row. The velocity contours are drawn in steps of 20 km s"' . 



saturated foreground stars in both images. The mean value of 
the flux ratio of on-line and the off-line filters yielded an initial 
scale factor of 0.053 with a standard deviation of 0.004. Then we 
created a set of continuum subtracted images with scale factors 
ranging from 0.045 to 0.061. A careful inspection of the images 
showed that the best scaling factor was between 0.049 and 0.053. 
We finally adopted a factor of 0.051 and estimate that it has an 
uncertainty of 7%. 



Observations of spectrophotometric standards were not 
available for our data. Our continuum subtracted image was 
therefore flux calibrated against previously publ ished Ha pho- 
tom etry for the H ii regions of this galaxy (iHodgd 119761: 
Ken nicutt & Hodge 1980). We i dentified a tot al of 26 Hn re- 
gions in both the finding charts of lHodgd d 19761) and our Ha im- 
age. The Hii regions were selected to be well distributed across 
the face of the galaxy and to cover a wide range in fluxes (from 
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^ 2 X 10"'"^ ergscm"^ to ^ 3.8 x 10"'"' ergscm"^). This com- 
parison yielded a calibration factor of (1.22 + 0.05) x 10"'^ 
ergs"' cm "^ for our Ha image, which assuming a distance of 
7.3 Mpc jSharina et al.l [T996I) . implies a luminosity of (7.8 + 
0.3) X lO^'* ergs"' per instrumental count. 

The Hi: region catalogue of NGC 628, i.e., the determina- 
tion of the position, size and luminosity of the H ii regions of 
the galaxy, was produced with the REGION so ftware package de- 
veloped by C. Heller (see IZurita et al]|200ll for an application 
of this code). The proce dure has been extensively described i n 
a number of papers (e.g.. lRozas et al.iri999l;lRelaiio et al.ll2005l) . 
This software has the advantage that it allows for the definition 
of local background regions above which the H ii regions are 
detected, and for manual editing of the preliminary results, al- 
lowing for changes if necessary. The selection criterion for con- 
sidering a feature in the Ha continuum subtracted image as an 
Hii region is that the feature must contain at least 16 pixels (an 
area approximately equal to the spatial resolution), each of one 
having an intensity of at least three times the r.m.s. noise of the 
image above the local background intensity level. This selection 
criterion defines the lower luminosity limit for the detection of 
an H II region in the image, which in this case corresponds to 
log Lho ~ 36.9 dex. With this selection criteria we detected a 
total of 376 Hii regions in NGC 628, with luminosities ranging 
from log Lhu - 36.9 to 39.3 and diameters up to 316 pc. Our 
measured fluxes agr ee within 10 - 20% with those reported by 
iKennicutt & Hodgd (Il980l) . An over plot of the individual H ii 
regions on our observed data is shown in Fig. |2] This catalogue 
has been made available on the CDS online service^ The CCD 
Ha image observed at La Palma is not shown in this article. 

4. Kinematic Analysis 

4.1. The Kinematic Analysis Principle 

The first step of our analysis is based on the assumption that 
the observed gas lies in an infinitesimally thin, and predomi- 
nantly rotating, circular disc. Observationally, discs show an ex- 
ponential light distribution, for which the equation for the rota- 
tiona l line-of-sight velocity component was derived by Freeman 
(I1970.) . When axisymmetric radial and vertical velocities are 
present, and taking into account the effect of projection and the 
convention that positive line-of-sight velocities correspond to re- 
cession, the line-of-sight velocities Vios can be written in polar 
coordinates as: 



V\os(R, lA, = Vsys + yrot(^) COS i// sin -H 

Vrad(^) sin tfr sin / + V-{R) cos /, 



(1) 



where Vsys is the systemic velocity of the galaxy, Vrot and Vr^d 
are the rotational and radial velocities, V, is the vertical velocity 
component, and / is the inclination at which the disc is projected 
on the sky. In the case of more inclined systems, where the im- 
pact of the radial and rotational velocities on Vios is dominant, 
the la st term of equation ([TJ can be ignored (e.g., Zuritaet al. 
12004 . However, in face-on galaxies and in cases where out of 
plane velocities are large, the effect of V, can in principle be 
measured. 

Gas kinematics in real galaxies exhibits radial and/or ver- 
tical motions due to, e.g., the presence of bars or spiral 



arms, which create angle-dependent velocities which cannot 
be described using equation ([T). Hence, additional ingredients 
in the analysis method ar e required. Franx et al. (1994) and 
ISchoenmakers et aP (1 19971) demonstrated that an observed ve- 
locity field can be decomposed into a harmonic series, the co- 
efficients of which provide quantitative information about non- 
axisymmetric structures: 



yios(^, (p, i) = Co 



^ c„cos(n 



0) + s„ sin(n0) sin i 



(2) 



^ The Hii region catalogue is made available in electronic format 
at the CDS via anonymous ftp to cdsarc.u-strasbg.fr or via 



[http : //cdsweb . u- strasbg . f r/ , 



where cq is the systemic velocity, ci is the circular velocity, 
and si corresponds to the axisymmetric radial velocity. The an- 
gle <p is measured in the plane of the orbit and is zero at the 
line of nod e s. Th is formalism has b e en successfully u sed b y 
'Wong e t al.l (|2004: iFathi et al.l (12 0051): 'E mseUem et akl (l2006l): 
Krainovi det al.l (l2006l) : lAnard et al. (2006) to model the effects 
of gravitational perturbations by bars and spiral arms on an ax- 
isymmetric velocity field. 



4.2. Kinematic Analysis Results 
4.2.1. The Velocity Field 

We started our analysis of the Ha velocity field by fitting an ex- 
ponential disc model to the observed velocity field. Fixing the 
centre to the kinematic centre position (see further down in this 
section), the best-fitting model yields an exponential disc scale 
length of 64±2", kinematic PA= 15+5°, and ysys=683+3 km s"' 
(see middle panel of Fig. [TJ. Although this fitting method also 
delivers the central surface density, which can be used to study 
the mass-to-light ratio of the disc, we find this topic to be beyond 
the scope of the current paper We investigated the robustness of 
this model by omitting regions where we expect large kinemati- 
cal discrepancies, such as the innermost regions or the outer parts 
of the arms where the residual velocities are large. We found that 
the fits delivered very similar parameters every time. Moreover, 
it is noteworthy that this procedure is computationally very ef- 
ficient as it delivers the kinematic parameters for all our 21952 
observed data points in 1 .2 seconds. A computational time much 
shorter than the analysis method described next. 

An alternative to fitting the entire field simultaneously 
is to section the field in concentric rings and to deri ve the 
kinematic parameters in a non-parametr ic fashion (Rotsi [T975t 
iBeauvais & Bofliunll200U lTeubenll2002h . We derived the kine- 
matics according to this method, for a simple rotating model, 
as well as the Vios according to the harmonic decomposition in 
equ ation dUl. The details of the fitting procedure are presented 
in Fathi et al.l (2005). The resulting best-fit velocity models are 
shown in the bottom panels of Fig.[T] In an attempt to fit the ve- 
locity field using equation ([T]i, we found that the radial velocity 
component is identical to the S[ term of equation whereas at 
all radii, V, stays at zero to within the errors. We conclude that 
the vertical velocities could not be determined from our data, 
either because they are too small to be detected with our obser- 
vations, or because they are mostly driven by small-scale pro- 
cesses, and thus the effects are smeared out within each ring. 

It has long been known that erroneous values for the posi- 
tion of the centre, the PA, or the inclination can cause recognis- 
able signatures in the residual velocity fields. Here we w ill not 
elaborate on this, but refer the reader to 'F athi et al.l (12005 ) for a 
detailed discussion. In the current study, we have made sure that 
none of the known systematic features appear in the residual ve- 
locity fields presented in Fig.[T] 
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Fig. 2. The Ha velocity dispersion (left) with over-plotted individual H n regions as described in section [373] The size of each circle 
indicates the size of the corresponding H ii region. The Ha intensity map with circles indicating individual H ii regions (right), and 
rectangles indicating approximate location of regions used to derive the diffuse interstellar gas properties (see section|5l). Spectra for 
one of these regions, obtained by co-adding the signal from 373 1 individual pixels, is illustrated. The centre of the galaxy is marked 
in both panels. 



Our derivation of the kinematic parameters is done in an iter- 
ative way. We first fitted the kinematic centre by calculating the 
average value of the (x, y) position of the centre. We then cal- 
culated the average PA, followed by the average Vjys, which we 
derived to be 19 + 7° and 682 ± 2 km s"' , respectively. The aver- 
ages were calculated usin g Tukey's bi-weight mean formalism 
jMosteller & Tukevlll977h which is particularly advantageous 
for being insensitive to outliers. Finally, we fitted the harmonic 
terms of equation (|2]). All these parameters and their profiles as 
a function of galactocentric radius are presented in Fig. |3] In 
NGC 628, one clear and one disturbed spiral arm and an oval 
distortion are the most prominent morphological features, the 
signatures of which we explore with our kinematical study. As 
pointed out by Schoenmakers et al. ( 1997), the first three terms 
of equation Q should deliver sufficient information for studying 
gravitational perturbation effects up to and including m - 2. 

A patchy gas distribution causes a non-uniform distribution 
of pixels, and it is important to make sure that all the rings in- 
clude enough pixels to fit harmonic parameters to a reasonable 
accuracy, especially in the innermost and outermost regions. We 
adopt a geometric increase of the ring radii with the radial width 
of the rings increasing by a factor 1 + step, where we adopt a 
step of 0.15. This corresponds to an innermost ring radius of 
6'.'5, containing 35 pixels. We only fit the inner 250" of the ob- 
served velocity field, since outside this radius the rings are not 
sufficiently covered and hence do not deliver reliably derived 
parameters. We found that every ring contained at least 33 data 
points, which normally leads to reliable results when fitting 6 
parameters simultaneously. 

We calculate the errors for the tilted-ring and harmonic pa- 
rameters indicated in Fig. [3] by means of Monte Carlo simu- 
lations. Repeated application of the tilted-ring method to the 
Gaussian-randomised gas velocity field yields the uncertainties 
on the harmonic parameters. Our 200 simulations show that 



choosing geometrically increasing ring radii with steps of 0.15 
indeed yields satisfactory errors. The errors are low, and the 
parameter profiles are represented by an adequate number of 
points. 



4.2.2. The Velocity Dispersion 

The large-scale velocity dispersion is one of the key parameters 
in the study of the physical state of the interstellar medium. For 
NGC 628, our selection criteria for trusting the Gaussian fits to 
the spectra allow us to measure reliably the gaseous cr intrinsic to 
the galaxy (see Section [T2] i. In the same fashion as carrying out 
the tilted-ring decomposition, we section the cr map in order to 
derive the average of this quantity as a function of galactocentric 
radius. Fig. |4] illustrates the average cr values in each galacto- 
centric radius bin, where we find that cr is almost constant at all 
radii with an average value of 16.7 + 1.8 km s"'. The average 
cr values reach a peak at around 1 10", followed by a slight de- 
crease outside this radius. The decrease is, however, well within 
the errors. 

The Hii region catalogue described in section [331 provides 
positions and sizes for individual H ii regions. In order to assess 
the total contribution of the H ii regions to the cr profile, we re- 
move all those pixels corresponding to H ii regions, and average 
over radial bins again. We calculate that excluding all H n re- 
gions, this quantity remains almost unchanged at cr = 16.4 + 1 .6 
km s The cr for the individual Hii regions also exhibits no 
trends as a function of galactocentric radius (see Fig.|4]i. We find 
that almost all radii within the galaxy include large as well as 
small H II regions, with an indication of a deficit of H ii regions 
at radii smaller than ^ 40", i.e., 1.9 kpc. 
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Fig. 3. The kinematic parameters derived from the velocity field of NGC 628, assuming an inclination of 6.5°. Radial profiles of 
the kinematic centre, PA, and systemic velocity are presented, the bi-weight mean averages (indicated by the horizontal dashed 
lines; mean values labelled when appropriate) of which have been calculated to constrain the geometry of the large-scale gaseous 
disc. The rotation curve and higher harmonic parameters have been derived for the obtained disc geometry, and the error bars were 
calculated as described in Section |4~2l For a brief discussion of the higher-order parameters (c2, S2, C3, 53), see Section|4] The bottom 
right panel shows the versus si diagnostics for the region inside the corotation radius of the spiral arms (open symbols), to which 
we have fitted a line. The slope of the line and its uncertainty is given at the bottom of the panel, where the negative slope indicates 
the effect of an m = 2 perturbation inside the corotation radius (see Section|53J. 



5. Discussion 

5.1. The Global Structure of NGC 628 

We must consider a number of global properties of NGC 628 be- 
fore building an evolutionary scenario for late-type spiral galax- 
ies of this type. 



It is well established that NGC 628 cannot have under- 
gone any e ncounter wit h satellites or other galaxies in the past 
10 years (IWakker & van Woerden 1991; Kamphuis & Briggs| 
1 1992 *). The presence of a large-scale elongate d structure, an oval , 
is confirmed from the isophotal analysis by iLaine et akl (l2002h 
at radii around 50" - 2.3 kpc. These authors also point out 
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InUlnsk: galaclocenlric radius ^arc^&c] 

Fig. 4. The velocity dispersion cr corrected for the instrumental 
dispersion as a function of radius (top panel), and the total cr for 
individual Hn regions (bottom panel). The vertical shaded area 
indicates the position of the nuclear ring-like structure. In each 
panel, the horizontal dashed line outlines the average value. 



signatures of a nuclear bar at nearly 2". This galaxy has one 
very prominent spiral arm to the south, with one or several dis- 
turbed spiral arms to the north, both extending from the nuclear 
regions out to well beyond the oval. The spiral filaments are pos- 
sibly disturbed by interaction with t he two large high-veloc ity 
clouds on either side o f the disc (e.g.,| Kamphuis & Briggs 1992; 
iBeckman et"a!1l2003h . ICepa & Beckma n (1990) found 'that the 
star formation efficiency is at its minimum at ^ 110", which 
they interpreted as the corotation radius for the spiral arms. In 
this galaxy, the O/H ratio has been found to decrease slightly 
from the centre outwards (Talent 1983|) , and the Ha luminos - 
ity decreases monotonically out to 450" dLelievre & Rovll2000h . 
The oxygen abundance decrease in NGC 628 is rather small, 
thus moderate mixing of the disc material driven by the oval dis- 
tortion is allowed (Zaritskv et al. 1994). 

Gas and stars are often found to be misaligned and/or dis- 
placed in spiral galaxies, and the detection of such misalign- 
ment has been eq ually successful in low-inclination and hi ghly 
inclined systems ( Fathil l2004t iFalcon-Barroso et"aL 1 120061) . In 
the case of an external origin of the gas, the nature of the 
orientation of the gas is thought to depend mostly on the 
initial angular moment um and mass of the accreted material 
jOuinn & BinnevI |1992|) . In this case, an observed kinematic 
misalignment can be interpreted as due to accreted gas rotating 
in a disc which is warped or incline d with respect to the stellar 
disc jSteiman-Cameron et alJlT992h . Analysing high resolution 
archival Hubble Space Telescope (HST), images in the optical 
and near infra-red, we have found that, within the uncertainties, 
the photometric PA and centre of NGC 628 are in agreement, to 
within six degrees and four arcseconds respectively, with the Ha 
kinematic values for these parameters (see Fig.©. 

Figure. |6] shows that the Ha rotation curve rises out to a 
galactocentric radius of 12 kpc, confir r ning p revious neutral hy- 
drogen observations by IWevers et alj (Il986l) (see also [T]i. The 
curve is quite bumpy, but shows a very clear velocity drop close 



Table 1. The derived deprojected rotation velocities and corre- 
sponding errors according to our tilted-ring decomposition. 



Radius [arcsec] 


Ko, [kms-'] 


Vrot error [km s ' ] 


6.5 


140.5 


40.5 


8.9 


115.0 


19.0 


10.5 


88.8 


23.3 


12.3 


81.7 


24.4 


14.4 


70.3 


19.9 


16.8 


61.6 


17.2 


19.7 


56.6 


17.3 


23.0 


79.9 


16.2 


27.0 


98.8 


11.6 


31.5 


117.6 


10.4 


36.9 


145.4 


9.2 


43.2 


128.7 


9.5 


50.6 


126.8 


7.3 


59.2 


131.9 


7.4 


69.2 


150.4 


7.3 


81.0 


156.5 


6.9 


94.8 


158.0 


5.7 


110.9 


155.6 


3.1 


129.8 


172.9 


4.2 


151.8 


140.5 


3.1 


177.7 


172.4 


5.1 


207.9 


177.3 


6.0 


243.2 


182.4 


3.4 



to around one kpc from the centre. This fall in the rotational ve- 
locity coincides quite closely radially with a dip in the B - V 
index. Although the non-parametric curve (dots with error bars) 
is not very well fitted by the exponential profile described in 
Section m we find that an exponential profile could reasonably 
describe the unperturbed disc component. The exponential mass 
distribution derived from our velocity field is almost exactly that 
derived from /?-band imaging bv .Natah et alj (Il992h . Our veloc- 
ity field yield s a scale length of ~ 64" which is 2.8 kpc on the 
galaxy, while iNatali et"aLl d 19921) used their photometric image 
to derive a scale length of 3.05 kpc, when their distance is con- 
verted to ours. 

Comparing only the kinematics of a non-parametric rotation 
curve with the exponential disc model is not straightforward for 
the case of NGC 628. In Fig. [1] we find that the residuals from 
the exponential disc model do not differ significantly from those 
yielded by the best fit rotating model. Both residual maps ex- 
hibit significant non-circular velocities, in particular in the re- 
gions surrounding the prominent spiral arm to the south, and 
around 50" to the north-east. Although we find that these resid- 
uals are taken out by fitting a harmonic decomposition (see bot- 
tom row in Fig.lTJ, the low inclination of the system, and hence 
the low ampUtudes of the non-circular motions, makes any de- 
tailed i nterpretation diffic ult. However, following the diagnos- 
tics by IWong et al] (|2004|) . our velocity field confirms that the 
general overall trend of the third versus the first sinusoidal term 
follows a negative slope inside the corotation radius (see bottom 
right panel of Fig. [3]), which suggests that an m = 2 perturbation 
is the most prominent in the observed velocity field. This nega- 
tive slope could be the kinematic confirmation of the presence of 
the oval perturbation. For t his diagnosis, we have u sed the coro- 
tation radius as derived by jCepa & BeckmanI 19901) . and assum- 
ing that the pattern speed of the spiral structure is the same as 
that of the oval structure, based on the response of galaxy disc to 
an analytic bar or oval potential (e.g., ISanders & Huntlevll 19761) 
as well as observational and theoretical results for galaxies such 
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Fig. 5. Zooming into the central regions of the velocity field of 
NGC 628. The HST //-band image contours (two innermost) 
and the FaNTomM continuum image contours (five outermost) are 
outlined. The kinematic centre and PA=19° are marked with a 
cross and the straight line, respectively. In this figure, Unear in- 
terpolation has been applied to fill the missing pixels for a clearer 
illustration. It is clear that the velocities drop at a radius around 
20", i.e., around the nuclear ring outlined by a dashed ellipse. 
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Fig. 6. Rotation curves assuming two different disc inclinations 
(top panel) together with an outline of the B - V colour profile 
from iNatali et al.l (Il992h (bottom panel). The solid curve indi- 
cates the rotation curve corresponding to a model of an expo- 
nential disc profile with a best fit scale length ~ 3 kpc, using 
6.5° inclination. See section |5T2l for further discussion. 



as NG C 1068 (lEmsellem et aHIIOOeh ESQ 507-16 (iBvrd et al.1 
[19981), and NGC 1365 (i Lindblad et aniT996h . 

The cr map presented in Fig. |2] displays local variations 
throughout the observed field, in comparison with the location 
of the individual Hii regions. We calculate the a of individ- 
ual H n regions by quadratically adding all the spectra within 
each H ii region, to which we then fit a single Gaussian profile to 
derive the kinematics. Analysing the dependence of the cr with 
galactocentric radius, shows that this quantity stays almost un- 
changed all the way out to 250", regardless of whether we take 
into account the individual Hii regions or not (see Fig.|4]i. This 
is not unexpected given that less than 18% of all pixels belong 
to individual Hii regions, and their average cr is not consider- 
ably higher than the mean cr throughout the disc, which is 16.7 
km s"'. 

The constancy with radius of the ionized Ha emit- 
ting gas veloc ity dispersion is ana logous to the constancy 
of cr for CO ([Combes & Becquaerll ,1997 ) and also for H i 
(iKamphuis & Briggslll992ir While the densities of the gaseous 
components are decreasing with radius, their dispersions vary 
very little. This is an interesting kinematic constraint on overall 
structure. The fact that the cr for the diffuse ionized component 
is larger than those of CO and H i is consistent with a larger scale 
height for this component, and this is as might be expected from 
detailed work in the Milky Way by Reynolds and his group (e.g., 
Madsen & Reynolds 2005). We will discuss this further in sec- 
tion |531 

We derived the cr of the diffuse interstellar by carefully se- 
lecting regions within the field with no significant Ho- emission 
line, i.e., with emission line amplitude-over-noise less than 3 
(see Fig. |2]l. The spectra from sections of these regions were 
then co-added quadratically, and a single Gaussian was fitted to 
the result. We confirm that the velocities for these selected re- 
gion are in agreement with Fig. [1] The emission line cr, on the 
other hand, was found to be similar to that of the ionized gas. 
For the diffuse gas, we derive the weighted average velocity dis- 
persion cr = 16.3 +1.3 km s This val ue is significantly lower 
than that found by many studies such as lWang et al.l (11997 ') who 
derived for five late-type spirals, the diffuse Ha average velocity 
dispersion larger than ^ 30 km s ' . The value we have derived 
shows that our data cannot be used to properly probe the diffuse 
interstellar medium. 



5.2. The innermost kiloparsec 

In Fig. |6l we see that the rotation curve of NGC 628 starts at 
high values, drops to a minimum just before a radius of 20", 
corresponding to 0.9 kpc (n.b. 1" ^ 47 pc), then rises again. 
Under normal conditions and mass-to-light ratios, such a central 
rise and fall is indicative of a rapidly rotating compact central 
structure, in this case smaller than 20", which fi ts well with the 
radius of the proposed nuclear ring (IWakker & A dler 1995). 

We recognise that the derivation of the kinematic centre, us- 
ing our large field of view, implies relatively large errors, of 8.6" 
and 9.2" in the x and y directions, respectively. Allowing for the 
kinematic centre to vary, the prominence of the circumnuclear 
rotation velocity rise changes slightly. We tested the rotation 
curve drop at around 15-20" by displacing the kinematic centre 
by up to 5 pixels (corresponding to 8"or 376 pc), after which we 
derived the rotation curve. For 10 different tests, we found that 
the location of the drop stays in the radius range 15-20", and 
that the amplitude of the drop-and-rise changes by less than 30 
km s ' (compare with Fig. |6] where the rotation curve drops by 
about 90 km s '). Thus, changing the centre to within the errors. 
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and re-calculating the rotation curve using the tilted-ring method 
yields consistent results which confirm that the innermost veloc- 
ity drop of the rotation curve is genuine. 

Centrally rotating structures have been observed in galax- 
ies, and are beUeved to form as a result of mass concentration 
caused by a large-scale perturbation or infalling satellites (e.g.. 



42.5 85.1 



Radius larcsec] 

127.6 170.2 212.7 255.3 



Whitehurst & Roberts 1972; Goad 1976; Bertola & Capaccioli 
1977 ; Shlosman&_Nog uchi 1993; Seifert & Scorza 1996). In 
the case of NGC 628, the clear isolation of this galaxy (e.g., 
iKamphuis & Briggslll992h complicates any external formation 
mechanism for the observed central mass concentration. The 
large-scale oval structure on the other hand, could transfer matter 
from the outer zones to regions nearer the centre of the galaxy. 
Alternatively, the inner disc-like component may be an observa- 
tional effect caused by a change of inclination between the inner 
and o uter component, i.e., an inner warp (e.g., Shen & SellwoodI 
12006 ). Thus the centre could be at higher inclinations, displaying 
higher deprojected line-of-sight velocities than if it was at the 
same inclination as the outer parts. As demonstrated in Fig. |6] 
the first velocity points of the rotation curve could be at a value 
lower than (or equal to) the R-2Q" velocity if the inner 20", cor- 
responding to 900 pc, w ould be at an inclination ^15 degrees 
larger than the outer disc. IWakker& Ad"ieil(ll995 h suggested that 
the nuclear molecular ring might be in a different plane than the 
outer disc. The nuclear CO deficiency detected by these authors 
could be only a relative deficiency, since only 45% of the single- 
dish flux was recovered in their interferom etric observat i ons. I n 
the higher spatial resolution observations bv lGanda et alj (l2006h . 
the nuclear ring appears much more circular, arguing in favour 
of the nuclear ring residing in the same plane as the outer disc. 
Moreover, the ellipticity profile of the high resolution HS T op- 
tical and near-infrared images argue against the presence of a 
nuclear warp in this galaxy. 

The evidence for the presence of a central structure in 
NGC 628 is not restricted to our Ha velocity fi el d. The position- 
velocity diagram presented in iDaigle et al.l (l2006al) clearly 
shows the presence of a rapidly rotating component in the cen- 
tre of this galaxy. In addition, our isophotal analysis of archival 
/-band and //-band HS T data confirm the presence of a nuclear 
component within the central 20" with a surface brightness pro- 
file distinctly deviating from that of the outer parts. Moreover, 
the B - V colour profile as well as the photometric study by 
iNatali etaTI (L1992) display the presence of an inner component 
(see our Fig. |6] as well as Fig. 7 of Natali et al. 1992). The 
latter authors interpreted the blue colour inside the ~ 20" ra- 
dius, as an indirect indication of infall of gas to the circumnu- 
clear regions. Investigating the archival images, w e also confirm 
the pr esence of the nuclear spiral arms found by Martini et al.l 
(I2003h . Moreover, as we describe in section 15.31 our kinematic 
data does indicate the presence of the structures necessary to in- 
duce infall of gas to the circumnuclear regions. 

We further investigated the SAURON gas and stellar kinemat- 
ics and distribution maps bv lGanda et al.l (12006'). and found con- 
firm ation for the 15 - 20" n uclear ring previously detected by, 
e.g.. lWakker & Adled (Il995h in their H/3 as well as [O m] emis- 
sion line distribution. The SAURON gas velocity field displays 
large line-of-sight velocities inside 15-20" radius as well as large 
up gas cr in the inner 15". The stellar component rotates in the 
same direction as the gas, and its velocity anti-correlates with the 
third harmonic term insi de this radius, indicating the signature 
of dis c-like kinematics ( lGerhardlll993l: Ivan der Marel & Franxl 
Il993h . The stellar cr displays a significant decrease of at least 30 
km s ' in the regi on inside 15", qualif ying NGC 628 as a cr-drop 
candidate (e.g.. lEmsellem et aLll2001l) . 




4 6 
Radius [kpc] 

Fig. 7. The angular frequency, Q, derived from our observed 
rotation curve (points with error bars), and that for a flat rota- 
tion curv e of 33 km (das hed line) corresponding to the curve 
drawn by lCepa & BeckmanI (1 1990) . Also drawn are the Q. + k/2 
and Q - k/2, the latter showing that the location of the nuclear 
ring is consistent with the location of the inner Lindblad reso- 
nance. The horizontal dotted line indicates the value of the pat- 
tern speed. 



A decrease of the stellar cr values in the centres of galax- 
ies requires a distinct dynamically cold and compact isother- 
mal nucleus in addition to the large-s cale stel l ar dis c, both with 
the same mass-to-light ratio (e.g., Bottemal Il989h . This pic- 
ture requires that substantial amounts of gas are directed to- 
wards the central regions, the gas cools and forms stars. The 
newly-born stars are thus concentrated toward the centre, and 
since they are formed from the low-dispersion gas compo- 
nent, their cr is lo wer than that of the o ld stellar component 
( Allard et al. 2005 | iBarbosa et alj 120061; iFalcon-Barroso et al.l 
2006c iFathi et al.l i2006). This is quite plausible under the cir- 
cumstances that gas accumulates in a nuclear disc or a bar, where 
dissipative cooling can be quite effective in producing the con- 
densation required to form stars. This ca n be a long-lived phe- 
nomenon jWozniak & Champavertll2006l) . 

The kinematic as well as photometric signatures observed in 
the innermost kpc of NGC 628 suggest the presence of a dy- 
namically cold disc-like structure which could be built by inflow 
from the outer parts towards the centre. 



5.3. Inflow, Circumnuclear Star Formation, and 
Pseudo-bulge Formation 

According to theory, the non-circular motions exerted by a bar 
or spiral arms funnel material towards the centre of their host 
galaxy. The gas accumulates in the inner Lindblad resonance 
region, where circumnuclear star formation can ensue, leading 
to a low luminosity starburst of up to about 10 Mg/yr (e.g., 
iKnapen et alJll995l ; lAthanassoulaii2005l and references therein). 

The corotation radius o f the s piral arms in NGC 628 was 
found by ICepa & BeckmanI (Il990h as the radius where the star 
formation efficiency is at its minimum. Using a projected flat ro- 
tation curve of 33 km s these authors derived a spiral pattern 
speed of 56 km s"' kpc"'. We investigate this value by deriving 
the angular frequency from our tilted-ring rotation curve (see 
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Fig. |7J. Considerable discrepancy between our observed curve 
and a flat one is found. For the corotation radius of 110" - 5.2 
kpc, our angular frequency curve corresponds to the density 
wave pattern speed D.p - 3 1 km s~' kpc~ , i.e . , alm ost a factor 
2 lower than the value from lCepa & BeckmanI (11990 ). An error 
of 15% in the presumed location of the corotation radius trans- 
lates to pattern speed errors of -1-5 and -3 km s ' kpc This 
is clearly due to our improved derivation of the rotatio n curve, 
different from that assumed bv lCepa & BeckmanI (Il990l) . 

We follow the observed signatures of an m = 2 perturbation 
and, using epicyclic approximation, derive the Q. + k/2 curves 
from our angular velocity curve. Given that this curve is rather 
bumpy, causing artifacts in calculating the numerical derivatives, 
we derive and present the smoothed Q. + k/2 curves. Figure [T] 
shows that our derived pattern speed together with the Q. - k/2 
implies the location of the inner Lindblad res onance radius at 
the p osition o f the nuclear ring as observed bv I Wakker & AdleJ 
(11995,) and bv lGanda et alj (12006) . However, the nuclear ring is 
not fully outlined in Ha, although two large Hii region com- 
plexes can be seen at its radius, the gas velocity dispersion rises 
slightly inside it, and there is a clear deficiency of H ii regions 
just outside it (see Fig.|4]i. 

Our Fabry-Perot observations of NGC 628 and the cur- 
rent analysis complement the already existing ingredients which 
allow us to build an evolutionary scenario for this galaxy. 
NGC 628 hosts a large-scale oval structure (or weak bar), one 
regular and a number of disturbed spiral arms, a nuclear ring 
at the inner Lindblad resonance radius, an inner, rapidl y rotat- 
ing, d isc-like structure, and a nuclear/se condary bar (Lai ne et alj 
2002), all signs of secular evolution dKormendv & KennicuttI 



20041) . This galaxy has, most likely, been exposed to some in- 



teraction more than 1 Gyr ago, which could cause the disturbed 
morphology in the south-eastern part of the galaxy. However, 
the time scale for this interaction is too short for its correspond- 
ing perturbation effects to reach the centre. The central parts of 
this galaxy display signatures of being built up by internal grav- 
itational effects of the prominent m - 2 perturbation, as a stel- 
lar cr drop is clearly obser ved. The y oung circumnuclear stel- 
lar populations observed bv ICornett et al. ( 1994), together with 
the dynamically cold centre invite the interpretation that we are 
observing the first phases of secular pseudo-bulge formation in 
NGC 628. 

Pseudo-bulges are components intermediate between 
spheroids and discs. They have shallower light profiles than 
spheroids, but steeper than discs. Galaxies with pseudo-bulges 
more frequently host starburs t s and circumnuclear star forma- 
tion dKormendv & KennicuttI l2004l) . Nuclear star formation 
in the form of nucle ar rings is also pre f erentially found in 
barred galaxies (e.g., iKnapen et al.ll2006l) . I Jo gee et al] d2005h 
have found that barred galaxies have large enough circum- 
nuclear molecular surface densities to trigger their observed 
circumnuclear starbursts. Molecular gas surface densities 
reach up to 20% of dynamical mass in the centre, eventually 
leading to the formation of a nuclear ring dPeeples & Martina 
|2006|) . Such concentration of gas and star formation, which 
are more common in late-type spiral galaxies, could lead to 
the formation of a pseudo-bulge (e.g., Fathi & Peletier 2003). 
The structural setup of NGC 628 fits well within this general 
picture, though this galaxy displays no complete circumnuclear 
star formation. Our photometric study of the archival data for 
NGC 628 suggests the presence of a central structure, which 
together with our kinematic analysis confirm that this central 
structure is disc-like and could be built up by inflow from large 
galactocentric distances. Similar cases, such as NGC 4593 and 



NGC 7690 dKormendv et alj2006l) . suggest that NGC 628 could 
be in the process of forming a pseudo-bulge. 

5.4. Velocity Dispersion for tlie Pliases of the Disc Gas 

Different gaseous components of the disc are expected to reach 
equilibrium. Different scenarios have been suggested for the for- 
mation and evolution of different disc layers ( Kroupal[2005L and 
references therein). The forces that govern the small scales are 
different from those that govern the large scale evolutionary pro- 
cesses. On the smaller scales, dissipation is expected to lower the 
gas cr, in order to reach dynamical equilibrium. The discrepant 
nature of the different gaseous components result in different 
cr levels, and thereby, different thicknesses. Observations have 
shown that the cr of the neutral molecular gas is ^ 7 km s"', cr 
of the warm neutral gas is ^ 14 km s"', cr < 30 km s"' for warm 
ionized Hcf-emit t ing gas, and cr > 30 km s ' for hot ionized gas 



(lRevnoldsl[ l98l 'KuUcarni & Fich''l98 5|). Hyd rodynamic simu- 
lations by de Avillez & Breitschwerdil 2005h have shown that 
these cr values remain constant for times as long as 0.5 Gyr, 
suggesting that thermal and turbulent pressures add to the to- 
tal pressure within the disc with constant rates. On large scales, 
gravitational instabilities form an essential heating mechanism, 
giving rise to vertical heating. 

Our analysis of the derived cr field of NGC 628 shines some 
new light on this subject. In Fig.|4] we find that cr reaches a peak 
value at arou nd 110", which cor r espon ds to the corotation radius 
suggested bvl Cepa & BeckmanI d 19901 ). followed by a slight de- 
crease outside this radius. Taking into account that the variations 
are within the errors, we find a near constancy of this param- 
eter with respect to galactocentric radius. We find that the av- 
erage overall cr does not vary when removing pixels belonging 
to individual H ii regions, clearly caused by the smearing effect 
of averaging along galactocentric rings. Nevertheless, our result 
confirms that the ionized gas cr is higher than its molecular CO 
and neutral counterparts. We interpret this result as the ionized 
gaseous component being thicker than the CO molecular and 
neutral gas components. The puzzling part of this result is the 
fact that removing the individual H ii regions does not change the 
average cr of the He-emitting gas. Could this be due to the highly 
ionizing effect of individual H ii regions on their surrounding in- 
terstellar medium? 

A substantial amount of star formation is expected in H ii re- 
gions, the global structure of which remains almost unchanged 
after a star-forming episode, implying the possibility to study 
the effect of feedback from stellar winds. Stellar winds not 
only increase the diffuse interstellar gas cr, bu t also contribut e 
to pushing it off the plane of the disc (e.g., Dettmar 2004|). 
One prominent observational example of this effect is the Milky 
Way, where the scale height of the diffu se interstellar ga s is as 
much as 10 times that of the stellar disc dRevnoIdsll 19891) . or in 
NGC 3938, where the ionised gas cr has been found to be sig- 
nifica ntly higher than that of the CO gas dJimenez- Vicente et al.l 
1 1999b . If these three galaxies are typical, they suggest that this 
effect could be common in spiral galaxies. However more de- 
tailed investigation is necessary, and in order to study the feed- 
back effect of the H ii regions on the diffuse interstellar medium 
in NGC 628, one needs to analyse the to tal emission-line pr ofiles 
within individual Hii regions (see, e.g.. lRelano et al.ir2005l) . 

It is crucial to quantify the vertical heating efficiency in 
NGC 628, along with the different emission-line ratios, which 
enables constraining the thickness of the different gaseous lay- 
ers. Furthermore, it is necessary to derive gradients in the ion- 
ization parameters along the sequence of transition between dif- 
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ferent states of the gas, and to investigate the underlying physi- 
cal changes in the ionized gas. Studies of this nature will allow 
quantifying the energy input from shock-excitation compared 
with that from photoionization. The currently available data set 
does not allow for analysis of this sort. 



6. Conclusions 

We have analysed Fabry-Perot interferometric observations of 
the Ha emitting ionized gas in the late-type spiral galaxy 
NGC 628 in order to study its internal kinematics over the full 
face of the galaxy with good angular resolution. Using a har- 
monic decomposition method we have detected the kinematic 
signature of an m = 2 perturbation and this, together with the 
form of the arms, and in broad agreement with previous kine- 
matic studies at lower resolution made with neutral gas, confirms 
that the galaxy is evolving with a minimum of external distur- 
bance. We have been able to derive the pattern speed, Q.p, of 
the spiral structure using the assumption that at corotation in the 
disc there is a minimum in the star formation efficiency, and find 
a value of 31^^ km s ' kpc This pattern speed value is also 
consistent with the location of the iimer Lindblad resonance ra- 
dius just below one kiloparsec, where a nuclear ring is observed, 
favouring the inflowing gas scenario according to theory. 

Using a continuum subtracted narrow band image in Ha we 
have identified 376 Hn regions with luminosities higher than 
36.9 ranging up to 39.3 dex, and with diameters up to 316 pc 
(see "Online Material" at the end). The lower limit here is ob- 
servational, conditioned by the angular resolution of the image 
and the distance of the galaxy. As star formation is very widely 
distributed within the disc we are not able to distinguish in the 
present exercise, at least for that part of the disc within the ra- 
dial range of the spiral arms, between true diffuse gas (as repre- 
sented in the Galaxy by the Reynolds layer) and a face-on "car- 
pet" of unresolved H ii regions. One clue here is the fact that 
the mean velocity dispersion for the ionized gas, 16.7 km s~*, 
(and almost constant over the radial range sampled, out to 12 
kpc), is hardly affected when the individually detected H ii re- 
gions are masked off, with a mean value in this case only 0.3 
km s"' smaller. This strongly suggests that in this face-on ob- 
ject with widely distributed star formation in the disc plane, the 
emission in Ha from the H n regions dominates any emission 
from the diffuse component. 

A key feature of the velocity field of NGC 628 is its regular- 
ity. The global effects of any asymmetries such as the oval dis- 
tortion detected are small, and in any case not very easy to detect 
in a face-on system. However we clearly detect kinematically an 
inner disc-like component in the inner kpc around the nucleus, 
a finding which is supported by corresponding variations in the 
B -V profile in this zone reported previously. There is at present 
no excess of star formation within this disc. Our kinematic anal- 
ysis suggest that the evolution of structure in NGC 628 has been 
driven by secular evolution of the disc in this galaxy. 
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Table 1. Complete list of the photometrically calibrated 376 individual H ii regions as described in section [331 The luminosities are given in units 
of erg s" ' , and the coordinates in Equ J2000. Note that here we only list a few of the brightest H ii regions. The entire catalogue is made available 
in electronic format at the CDS via anonymous ftp to cdsarc . u-strasbg . fr or via http://cdsweb . u-strasbg . fr/ 
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